NASA/TM— 20 1 0-2 1 6344 



Effect of High-Fidelity Ice Accretion Simulations on 
the Performance of a Full-Scale Airfoil Model 


Andy P. Broeren and Michael B. Bragg 
University of Illinois, Urbana, Illinois 

Harold E. Addv, Jr. 

Glenn Research Center, Cleveland, Ohio 
Sam Lee 

ASRC Aerospace Corporation, Cleveland, Ohio 

Frederic Moens 
ONERA, Meudon, France 

Didier Guffond 
ONERA, Chatillon, France 


June 2010 


NASA STI Program ... in Profile 


Since its founding, NASA has been dedicated to the 
advancement of aeronautics and space science. The 
NASA Scientific and Technical Information (STI) 
program plays a key part in helping NASA maintain 
this important role. 

The NASA STI Program operates under the auspices 
of the Agency Chief Information Officer. It collects, 
organizes, provides for archiving, and disseminates 
NASA’s STI. The NASA STI program provides access 
to the NASA Aeronautics and Space Database and 
its public interface, the NASA Technical Reports 
Server, thus providing one of the largest collections 
of aeronautical and space science STI in the world. 
Results are published in both non-NASA channels 
and by NASA in the NASA STI Report Series, which 
includes the following report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant phase 
of research that present the results of NASA 
programs and include extensive data or theoretical 
analysis. Includes compilations of significant 
scientific and technical data and information 
deemed to be of continuing reference value. 
NASA counterpart of peer-reviewed formal 
professional papers but has less stringent 
limitations on manuscript length and extent of 
graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary or 
of specialized interest, e.g., quick release 
reports, working papers, and bibliographies that 
contain minimal annotation. Does not contain 
extensive analysis. 

• CONTRACTOR REPORT Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. Collected 
papers from scientific and technical 
conferences, symposia, seminars, or other 
meetings sponsored or cosponsored by NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, often 
concerned with subjects having substantial 
public interest. 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific and 
technical material pertinent to NASA’s mission. 

Specialized services also include creating custom 

thesauri, building customized databases, organizing 

and publishing research results. 

For more information about the NASA STI 

program, see the following: 

• Access the NASA STI program home page at 
http://www.sti.nasa.gov 

• E-mail your question via the Internet to help@ 
sti.nasa.gov 

• Fax your question to the NASA STI Help Desk 
at 443-757-5803 

• Telephone the NASA STI Help Desk at 
443-757-5802 

• Write to: 

NASA Center for AeroSpace Information (CASI) 
7115 Standard Drive 
Hanover, MD 21076-1320 


NASA/TM— 20 1 0-2 1 6344 



Effect of High-Fidelity Ice Accretion Simulations on 
the Performance of a Full-Scale Airfoil Model 


Andy P. Broeren and Michael B. Bragg 
University of Illinois, Urbana, Illinois 

Harold E. Addv, Jr. 

Glenn Research Center, Cleveland, Ohio 
Sam Lee 

ASRC Aerospace Corporation, Cleveland, Ohio 

Frederic Moens 
ONERA, Meudon, France 

Didier Guffond 
ONERA, Chatillon, France 


National Aeronautics and 
Space Administration 


Glenn Research Center 
Cleveland, Ohio 44135 


June 2010 


Acknowledgments 


The success of this collaborative research project is due to the combined efforts of many individuals at NASA, ONERA and 
University of Illinois. In particular, the staff of the ONERA F 1 facility provided excellent support in the planning and execution 
of the aerodynamic testing. The ice simulations attached to the model during the testing were skillfully fabricated by the staff of 
the Aerospace Engineering Department Machine Shop at the University of Illinois. The authors at the University of Illinois were 
supported under NASA grant NCC3-1039 and ASRC Aerospace Corporation subcontract AS3304 from NASA Glenn Research 
Center. This support is gratefully acknowledged. 


Level of Review. This material has been technically reviewed by technical management. 


Available from 


NASA Center for Aerospace Information 
7115 Standard Drive 
Hanover, MD 21076-1320 


National Technical Information Service 
5301 Shawnee Road 
Alexandria, VA22312 


Available electronically at http://gltrs.grc.nasa.gov 


Effect of High-Fidelity Ice Accretion Simulations on the 
Performance of a Full-Scale Airfoil Model 


Andy P. Broeren and Michael B. Bragg 
University of Illinois 
Urbana, Illinois 61801 

Harold E. Addy, Jr. 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

Sam Lee 

ASRC Aerospace Corporation 
Cleveland, Ohio 44135 

Frederic Moens 
ONERA 

Meudon, France F92190 

Didier Guffond 
ONERA 

Chatillon, France F92322 


Abstract 

The simulation of ice accretion on a wing or other surface is 
often required for aerodynamic evaluation, particularly at 
small scale or low-Reynolds number. While there are 
commonly accepted practices for ice simulation, there are no 
established and validated guidelines. The purpose of this 
article is to report the results of an experimental study 
establishing a high-fidelity, full-scale, iced-airfoil 
aerodynamic performance database. This research was 
conducted as a part of a larger program with the goal of 
developing subscale aerodynamic simulation methods for iced 
airfoils. Airfoil performance testing was carried out at the 
ONERA FI pressurized wind tunnel using a 72-in. 
(1828. 8-mm) chord NACA 23012 airfoil over a Reynolds 
number range of 4.5xl0 6 to 16.0xl0 6 and a Mach number 
range of 0.10 to 0.28. The high-fidelity, ice-casting 
simulations had a significant impact on the aerodynamic 
performance. A spanwise-ridge ice shape resulted in a 
maximum lift coefficient of 0.56 compared to the clean value 
of 1.85 at Re = 15.9* 10 6 and M = 0.20. Two roughness and 
streamwise shapes yielded maximum lift values in the range of 
1.09 to 1.28, which was a relatively small variation compared 
to the differences in the ice geometry. The stalling 
characteristics of the two roughness and one streamwise ice 
simulation maintained the abrupt leading-edge stall type of the 
clean NACA 23012 airfoil, despite the significant decrease in 
maximum lift. Changes in Reynolds and Mach number over 
the large range tested had little effect on the iced-airfoil 
performance. 


Nomenclature 

a Airfoil angle of attack 

a sta u Stalling angle of attack, coincident with 

c Airfoil chord length 

Cj Drag coefficient 

C/ Lift coefficient 

C/ >a Lift-curve slope 

C/ max Maximum lift coefficient, coincident with a sta n 

C m Quarter-chord pitching-moment coefficient 

k Ice-roughness height or thickness 

L WC Liquid Water Content 

M Freestream Mach number 

MVD Median Volumetric Diameter 

Re Reynolds number based on chord 

x Chordwise position along airfoil 

y Normal position from airfoil chord line 

1.0 Introduction 

In many applications, simulation of the ice accretion 
geometry on a wing or other surface is required for 
aerodynamic evaluation. The measurement of aerodynamic 
performance from an aircraft or wing with actual ice 
accretions is very difficult and expensive and seldom a 
practical solution. Aircraft performance data with ice accreted 
in flight are available (Refs. 1 to 3), but limited, due to the 
cost and difficulty in acquiring these data. In flight, it is 
difficult to document the ice accretion geometry accurately 
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where shedding often occurs. The accretions are difficult to 
attribute to a particular icing condition due to the natural 
atmospheric variations through which the airplane traverses. 
For accretions formed in an icing wind tunnel the inability to 
provide adequate instrumentation, nonuniformity of the cloud 
over the model, poor aerodynamic flow quality, and cost also 
make aerodynamic measurements difficult and rare. The most 
common way to acquire iced-airfoil and wing data is to use 
simulated, artificial ice accretions in a dry-air wind tunnel, or 
in flight. This process requires methods for developing the ice 
accretion simulations. 

The best current technology for accurate ice accretion 
simulation is the mold and casting method developed at the 
NASA John H. Glenn Research Center (Ref. 4). In this 
method molds are made from an ice accretion generated in an 
icing wind tunnel, such as the NASA Glenn Icing Research 
Tunnel (IRT). From these molds, ice accretion castings are 
made that maintain the major features of the ice including the 
detailed surface roughness and the spanwise and chordwise 
variations. These castings are then instrumented to obtain 
high-fidelity aerodynamic data (Refs. 5 to 7) This is an 
expensive process and not practical in many situations, but 
does generate benchmark data for iced-airfoil and wing 
research. 

While much of our understanding of ice accretion 
aerodynamics is anchored by cast ice simulation data, there 
have been very few studies employing full-scale ice accretion 
cast shapes aerodynamically tested at full-scale Reynolds and 
Mach number. The second FAA/NASA Tailplane Icing 
Program provides full-scale, cast-ice performance data at full- 
scale Reynolds and Mach number, but the ice shapes were 
applied to a fully three-dimensional tailplane model (Ref. 8). 
The NASA Modern Airfoils Program (Ref. 9) currently 
provides the best airfoil performance data set at flight 
Reynolds and Mach numbers, but due to size limitations in the 
NASA Langley Low-Turbulence Pressure Tunnel (LTPT) 
used for the aerodynamic testing, the ice accretions were 
obtained on a 36-in. (914.4-mm) chord airfoil. Ice accretion 
scaling methods can provide scaling of the gross ice shape for 
moderate scales (Ref. 10), but not of the roughness and ice 
shape details known to be aerodynamically important in some 
cases. 

The simulation of ice accretions on airfoils and wings is 
needed or desirable for many applications. For flight or wind- 
tunnel tests of iced airfoils for certification or research, cast 
ice shapes are seldom available and lower-fidelity simulations 
are required. Computational methods must model ice on 
airfoils and wings and it is either impossible or impractical to 
model a full, rough 3-D accretion. Simpler geometries that 
accurately represent the key aerodynamic features are needed. 
Currently, there is only a limited understanding of how to 
accurately simulate the aerodynamic effect of ice accretion on 
lifting surfaces. Accurate simulations need to be based on an 
understanding of the fundamental flow physics for the 
simulation to be robust under a wide range of geometries and 


flow conditions. While there are commonly accepted practices 
supported by some research in this area, there are also 
remaining questions about the accuracy of such simulations 
(Ref. 11). 

Many of these open issues are directly addressed in the 
“Airfoil Ice-Accretion Aerodynamics Simulation” research 
program sponsored by NASA and ONERA (Ref. 12). The 
completed program will result in validated scale model 
simulation methods that produce the essential aerodynamic 
features of the full-scale iced airfoil. This program, 
summarized in Section 2.0 below, required the development of 
a full-scale, iced-airfoil aerodynamic database for validation 
of the subscale model simulation methods. These data also 
provide, for the first time, aerodynamic data for high-fidelity 
ice simulations on a full-scale airfoil model at flight Reynolds 
number. The castings were fabricated from ice accretion 
testing conducted in the NASA IRT and were applied to the 
leading edge of a 72-in. chord (1828. 8-mm) NACA 23012 
airfoil model. Airfoil performance testing was carried out at 
the ONERA FI pressurized wind tunnel over a Reynolds 
number range of 4.5xl0 6 to 16.0xl0 6 and a Mach number 
range of 0.10 to 0.28. A total of six different ice accretion 
casting simulations were tested. The purpose of this article is 
to present these results. 

2.0 Background: Iced-Airfoil 
Aerodynamic Simulation 

The full-scale, high-fidelity, iced-airfoil aerodynamic 
database reported in this article was motivated by the problem 
of aerodynamic simulation of icing effects. Bragg et al. 
(Ref. 12) describe this six-phase research program. The plan 
included both subscale and full-scale ice accretion and 
aerodynamic testing. By using subscale testing whenever 
possible to minimize the expense of full-scale testing, and as a 
way to validate the scaled ice accretion simulations, an 
efficient and cost-effective program was developed. In Phase 
1: Ice-Shape Classification, a review of the existing iced- 
airfoil aerodynamic literature was performed in order to 
classify ice shapes according to their aerodynamic effects. The 
objective of Phase 2: Subscale Model Ice Accretion Testing 
was to obtain high-fidelity ice shapes for the subscale model 
having the characteristics developed in Phase 1. In Phase 3: 
Subscale Model Aerodynamic Testing, the high-fidelity ice 
shapes obtained in Phase 2 were used to develop aerodynamic 
simulation methods on the subscale model. The objective of 
Phase 4: Full-Scale Model Ice Accretion Testing was to obtain 
high-fidelity ice shapes for the full-scale model having the 
characteristics developed in Phase 1. In Phase 5: Full-Scale 
Model Aerodynamic Testing , the ice-shapes acquired in Phase 
4 were used for aerodynamic testing on the full-scale model to 
obtain a benchmark data set for the validation of subscale 
simulation methods. Finally, in Phase 6: Simulation 

Validation Testing , the methods developed in Phase 3 were 
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used to scale and simulate the full-scale ice shapes for testing on 
the subscale model at lower Reynolds number. The objective of 
this phase was to “close the loop” by using the subscale model 
to reproduce the aerodynamic effects of the ice shapes tested on 
the full-scale model at high Reynolds number. Each phase of 
this research program is summarized below. 

2.1 Phase 1: Ice-Shape Classification 

The first phase of the research was to develop ice-shape 
classifications based upon the iced-airfoil aerodynamics. 
Bragg et al. (Refs. 13 and 14) defined four categories of ice 
accretion: roughness, horn, streamwise, and spanwise ridge. 
Roughness is defined as the type of ice accretion that affects 
the boundary-layer transition process and extracts momentum 
from the boundary layer. This leads to premature trailing-edge 
separation and increased drag and reduced lift. Specific 
examples of roughness and the other types of ice are discussed 
in Section 3.0. Horn ice is a larger accretion that has at least 
one protrusion orientated at a significant angle to the 
oncoming flow and is located in the airfoil leading-edge 
region such that the stagnation point is on the ice accretion. 
This results in a large separation bubble aft of the accretion 
that governs the aerodynamics. The separation bubble 
significantly increases the drag and reduces the stalling angle 
of attack and maximum lift. Streamwise ice is more conformal 
to the leading edge and thus does not result in a large 
separated flow region. There may be localized separation on 
the scale of the ice shape owing to discontinuities at the 
ice/airfoil junction. The boundary layer is significantly 
affected by the ice, resulting in earlier transition and trailing- 
edge separation. Spanwise-ridge ice is characterized by a large 
separation bubble aft of the ridge, but differs from horn ice in 
that ridges are located downstream of the leading edge. The 
airfoil upstream of the ice ridge is smooth such that the 
stagnation point is not located on the ice accretion, and a 
boundary layer develops upstream. This leads to an additional 
separation bubble upstream of the ridge thus making the three- 
dimensional characteristics of the ridge more important than in 
the horn ice case. It is important to note that there is often 
overlap among these classifications and that real ice accretion 
could be characterized by a combination of these descriptions. 
However, this type of flowfield-based understanding is useful 
from the simulation perspective. Thus, it was this 
understanding that guided the selection of full-scale ice 
accretions for the aerodynamic testing presented in this article. 

2.2 Phase 2: Subscale Model Ice Accretion 
Testing 

The objective of this phase was to obtain high-fidelity ice 
shapes having the characteristics of the four types described 
in Phase 1. To obtain these shapes, a subscale, 18-in. 
(457.2-mm) chord NACA 23012 airfoil model was subjected 
to simulated, in-flight icing conditions in the NASA IRT. The 


icing conditions were selected to produce accretions 
representative of the four classifications. Molds were made of 
selected ice accretions that were used to develop castings 
(Ref. 4) for aerodynamic testing in Phase 3. Blumenthal 
(Ref. 15) provides more details and results of this phase of the 
program. 

2.3 Phase 3: Subscale Model Aerodynamic 
Testing 

The objective of this phase was to develop simulation 
methods to reproduce the aerodynamics of the iced-airfoil. 
Each of the four classifications has different fundamental 
aerodynamics. If techniques can be developed to simulate 
these four shapes, then most ice accretions can be simulated. 
A variety of simulation methods already exist in the literature 
and in practice and these were categorized for subscale testing. 
The simulation categories considered were: 3-D castings, 2-D 
smooth and simple-geometry. To account for roughness 
associated with ice accretion, distributed grit-type roughness 
may be added to the non-casting simulation methods. The 
casting was considered to be the highest fidelity simulation. It 
is a benchmark for aerodynamic testing that the other lower- 
fidelity simulation methods are evaluated against. 

Numerous experimental trials have been conducted with the 
ice accretions obtained in Phase 2 to investigate the 
effectiveness of the various simulations and quantify the 
associated accuracy. These results are documented in several 
papers (Refs. 16 to 20) and are summarized here. The 
combined results indicate that different methods should be 
used to best represent the aerodynamics of the iced-airfoil 
according to the different classifications of ice accretion. In 
the horn-ice case, the 2-D smooth and simple-geometry 
simulations were nearly equal in effectiveness owing to the 
dominance of the separation bubble flowfield. The addition of 
roughness reduced the simulation fidelity when compared to 
the casting results. In contrast, the addition of roughness to the 
2-D simulations was required to obtain the best agreement in 
performance with the casting configuration for the streamwise 
ice accretion. In this case, the size and concentration had 
significant effects such that these parameters should be taken 
into account. The performance effects of the spanwise-ridge 
ice accretion was best reproduced by the 2-D smooth 
simulation. The spanwise-ridge shape did not result in a large, 
closed-two dimensional separation bubble and the resulting 
flowfield was highly three dimensional with evidence of 
trailing-edge separation. These effects complicated the ability 
to simulate the aerodynamics of the iced airfoil with lower- 
fidelity geometries. 

2.4 Phase 4: Full-Scale Model Ice Accretion 
Testing 

The objective of this phase of the program was to obtain 
high-fidelity ice shapes in each of the four categories on a full- 
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scale airfoil. Since this effort was designed to develop 
subscale simulation methods, it was important to have a set of 
benchmark data and firm knowledge of the full-scale 
phenomena. This information provides validation data for both 
the subscale testing and simulation methods developed. In 
order to obtain frill-scale aerodynamic effects information for 
the four classifications of ice, accurate representations of the 
ice were needed. Therefore, a full-scale, 72-in. (1828. 8-mm) 
chord, NACA 23012 airfoil, ice accretion model was built and 
tested in the NASA IRT. This effort provided the full-scale, 
reference ice shapes for the program. Molds were made of 
selected ice accretions and were used to develop castings 
(Ref. 4) that were used for the full-scale aerodynamic testing 
and are described in Section 3.0. 

2.5 Phase 5: Full-Scale Model Aerodynamic 
Testing 

In this phase, aerodynamic testing was conducted to 
establish a set of high-fidelity benchmark data useful for the 
development of simulation methods at smaller scale. Airfoil 
performance testing was performed in a pressurized wind 
tunnel. A total of six different ice accretion casting 
simulations were tested. As was the case for the subscale 
testing, ice shapes were selected according to the four 
classifications developed in Phase 1. In contrast to the 
subscale testing, only the high-fidelity ice casting simulations 
were tested on the full-scale model. This approach increased 
the number of different ice shapes that could be tested with the 
given amount of tunnel time. The purpose of this article is to 
present these results. 

2.6 Phase 6: Simulation Validation Testing 

In this final phase of the program, the ice shapes selected 
for the full-scale aerodynamic testing were scaled and 
simulated for testing on the subscale model. Lower-fidelity 
simulations of the full-scale ice shapes were developed and 
tested on the quarter-scale model at low-Reynolds number 
(1.8x1 0 6 ). Busch et al. (Ref. 21) quantified the accuracy of the 
subscale simulations against the full-scale data reported in this 
article. For horn ice, it was determined that accurately 
representing the height, angle, and location of the upper- 
surface horn using a 2-D smooth simulation provided an 
accurate representation of C/ max . Similar results were obtained 
for spanwise-ridge ice — the simulation of geometric details of 
the ice shape including roughness had only minor effects on 
the iced-airfoil C/ max . For roughness and streamwise ice, it was 
found that geometric scaling of the roughness height tended to 
cause unrepresentatively large penalties to Q,„ m and Q when 
the roughness was applied at high concentration. Reducing 
roughness concentration would decrease these penalties, but at 
this time no accurate method exists for measuring roughness 
concentration on an ice accretion. Therefore, modeling both 


roughness height and concentration is important in subscale 
iced-airfoil simulation for some types of ice. 

3.0 Experimental Methods 

All of the aerodynamic testing reported here for Phase 5 
was performed at the ONERA F 1 full-scale, pressurized wind- 
tunnel facility (Ref. 22). The closed-return wind tunnel has a 
test-section measuring 138-in. (3500-mm) high by 177-in. 
(4500-mm) wide by 433-in. (11000-mm) long. The maximum 
test section Mach number is 0.36 and the maximum stagnation 
pressure is 57 psia (3.85 bar). The unit Reynolds number can 
be varied up to a maximum of 6.0 X 10 6 /ft at Mach = 0.23. 
Total temperature is maintained via a heat exchanger located 
in the second diffuser downstream of the fan. The fan operates 
at constant speed while the test section Mach number is 
controlled by adjusting the pitch of the blades. The test-section 
inlet flow is conditioned through a 7.18-to-l contraction 
containing honeycomb flow straightener and three turbulence 
reduction screens. 

The 72-in. (1828. 8-mm) chord NACA 23012 airfoil model 
was mounted vertically in the test section as shown in 
Figure 1. The model span was 137.48-in. (3492-mm) and was 
mounted in the floor force balance. Small gaps between the 
bottom of the model and the test-section floor as well as the 
top of the model and the test-section ceiling were maintained 
so as not to cause mechanical hysteresis in the force-balance 
measurements. The model had a main chordwise row of 72 
pressure taps located at 43 percent span measured from the 
test-section floor (cf. Fig. 1). In addition, there was a row of 
20 taps oriented spanwise at x/c = 0.70 on the upper surface. 
The model was designed and built with full-span removable, 
interchangeable, leading-edge sections. The baseline leading 
edge had the clean NACA 23012 profile, while the alternate 
leading edge had a truncated nose geometry. The latter design 
facilitated mounting of the various ice-shape casting 
simulations. Accommodations were also made in the pressure 
tapping to allow for rapid connection of pressure 
instrumentation in the cast ice shapes. Also shown in Figure 1 
is the wake rake located one chord length downstream of the 
model trailing edge. The wake rake had 100 stagnation 
pressure probes spaced 0.79-in. (20-mm) apart and was 
located at a fixed spanwise station at 57 percent span above 
the test-section floor. 

Data acquisition runs were performed in angle of attack 
sweeps for increasing and then decreasing angle of attack at a 
constant sweep rate of 0.1°/sec. Data were also acquired at 
fixed angle of attack for selected angles over the range of the 
sweep and repeat runs were performed to ensure data 
precision. The data shown in this paper are for increasing 
angle of attack sweeps and have been averaged to the nearest 
0.5° in post-processing. The averaging method divided the 
data into 0.5° bins and averaged the results. For example, data 
from a = 2.75° to 3.25° were average to create a data point 
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x/c 


Figure 1. — Photograph of NACA 23012 airfoil model installed 
in ONERA FI wind-tunnel test section and model pressure 
tap layout. 

close to a = 3.0° The actual angle of attack value was the 
arithmetic mean of the angles of attack collected in that bin. 
During the sweeps, data were acquired from the force balance 
and tunnel conditions analog transducers. The model surface, 
test-section sidewall, and wake-rake pressures were acquired 
using an electronically scanned pressure system. The 
acquisition of these data were synchronized in time 
corresponding to the angle of attack sweep rate. Lift and 
pitching-moment coefficients were calculated from the force 
balance and from the integration of the surface pressure 
measurements. Good agreement between the integrated- 
pressure data and the force-balance data were obtained. In this 
paper, the lift and pitching moment data reported for the clean 
configuration were obtained from the surface pressures while 
the data reported for the iced configurations were obtained 
from the force balance. The force balance data are reported 
because, in many cases, the stall of the iced-airfoil 
configurations was characterized by unsteady flow. The 
available signal conditioning for the force-balance data 
allowed for more effective filtering of these unsteady effects. 
Drag coefficient was calculated from the wake pressures using 
standard momentum-deficit methods and these values are 
reported in this paper for all configurations. The performance 
coefficients were corrected for wind-tunnel wall effects using 
the methods of Allen and Vincenti (Ref. 23). The angle of 
attack sweeps were performed for a large range of Reynolds 
numbers and Mach numbers as shown in Table 1. The matrix 
was designed to isolate the independent effects of these 
parameters. Therefore, Reynolds numbers variations were 
performed at constant Mach numbers of 0.10 and 0.20 while a 
Mach number variation was performed at a constant nominal 
Reynolds number of 12xl0 6 . Due to operational constraints, 


the Reynolds number was not maintained entirely constant for 
each of these conditions, and therefore Table 1 indicates the 
variation over the course of the entire test campaign. 


TABLE 1.— MATRIX OF REYNOLDS AND 
MACH NUMBER CONDITIONS 


Reynolds 
number range 

Mach number 

0.10 

0.20 

0.28 

4.5xl0 6 to4.6xl0 6 

X 



8.0xl0 6 to 9.1xl0 6 

X 

X 


12.0xl0 6 to 12.3x10 s 

X 

X 

X 

15.8x10 s to 16.0x10 s 


X 



The experimental uncertainty in the performance 
coefficients was estimated using the methods of Kline and 
McClintock (Ref. 24) and Coleman and Steele (Ref. 25) for 
20:1 odds. Table 2 lists these uncertainties for both integrated- 
pressure and force-balance measurements, before the wall 
corrections were applied. The values were calculated based 
upon the clean model configuration at Re = 8.1xl0 6 and M = 
0.20. The uncertainties are expected to be identical for the 
iced-model configurations. The absolute uncertainties in 
Table 2 are inversely proportional to the dynamic pressure 
(except for a). This condition was selected because it 
corresponds to the average dynamic pressure over the range of 
conditions (Table 1). Therefore, conditions having lower 
dynamic pressure would have slightly larger uncertainties 
while conditions with higher dynamic pressure would have 
slightly lower uncertainties. All of these uncertainties were 
acceptable for the purposes of this investigation. The relative 
uncertainty in C m (both pressure and balance) seems large for 
this example because of the small reference value. For cases 
where the C,„ values were larger, e.g. in the iced-airfoil case, 
the absolute uncertainty would be similar, therefore resulting 
in a lower relative uncertainty. This is also the case for the 
uncertainty in drag coefficient. Several repeat runs were 
performed for both clean and iced configurations and the run- 
to-run variations in the coefficients were much smaller than 
the uncertainties listed in Table 2. The good agreement 
between the integrated-pressure and force-balance data 
provided further assurances of data quality. 


TABLE 2.— ESTIMATED EXPERIMENTAL UNCERTAINTIES. 


Aerodynamic 

quantity 

Reference 

value 

Absolute 

uncertainty 

Relative 

uncertainty 

a 

8.01° 

±0.02° 

±0.25% 

Ci Balance 

1.095 

±0.010 

±0.93% 

C,„ Balance 

-0.0144 

±0.00071 

±4.9% 

c P 

-1.057 

±0.015 

±1.41% 

Ci Pressure 

1.096 

±0.0070 

±0.64% 

C m Pressure 

-0.0148 

±0.0024 

±16.5% 

Q Wake 

0.0086 

±0.00048 

±5.5% 
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A total of six ice accretion simulations were tested. Using 
the ice accretion classifications developed by Bragg et al. 
(Refs. 13 and 14) there was one horn shape, one spanwise- 
ridge shape, two roughness shapes and two streamwise shapes. 
The casting simulations were fabricated from ice accretion 
moldings acquired during Phase 4 of the program as described 
in Section 2.0 and in more detail by Bragg et al. (Ref. 12) 
Table 3 summarizes the icing conditions associated with the 
ice accretions selected for aerodynamic testing. These ice 
accretions are documented in Figures 2 to 7. The horn shape in 
Figure 2 shows the classic upper-surface horn typical of this 
glaze -type accretion. The streamwise shape in Figure 3 is 
more conformal to the leading-edge radius with a smooth zone 
on the nose followed by downstream rime feather roughness. 
The glaze roughness in Figure 4 had a smooth zone in the 
stagnation region followed by large roughness downstream. 
The spanwise-ridge shape in Figure 5 was formed by applying 
a heater mat to the model leading edge. The heat input and 
icing conditions were adjusted to accrete the upper and lower 
surface ridges shown in Figure 5. The streamwise shape in 
Figure 6 was selected for aerodynamic testing because it had a 
more “pointed” geometry at the leading edge and was less 
conformal to the surface than the other streamwise shape 
(Fig 3). The remaining roughness case in Figure 7 was formed 
in cold conditions resulting in very fine rime feathers that 
were very different in size and distribution from the glaze 
roughness case in Figure 4. 

For each ice accretion a set of casting simulations were 
fabricated for aerodynamic testing. The castings were made 
from a hard two-part, polyurethane-type material of very low 
viscosity required to flow into the intricate details of the ice 
molds. The moldings of the ice accretions were limited to a 
span of 15-in. (381-mm) at the center of the IRT model. This 
was sufficient length to reproduce the spanwise variation in all 
of the ice accretion geometries. Therefore, multiple sections of 
casting were fabricated for each ice shape to cover the 137.48- 
in. (3492-mm) span of the NACA 23012 model. A total of 10 
sections were used, each having a finished length of 13.75-in. 
(349.3-mm). The sections were bolted on to the alternate 
leading edge of the model and the interfaces were sealed with 
silicone. A completed installation is shown in Figure 8. One 
casting section was instrumented with static pressure taps. 
This was located at the same spanwise station as the main 
chordwise row on the model (43 percent span measured 


from the test-section floor). The orifices were drilled directly 
into the casting and located in areas conducive to good surface 
static pressure measurements based upon previous research 
(Ref. 15). The number of pressure taps on the entire 
removable section ranged from 29 to 36 with some ice shapes 
requiring higher resolution. The pressure orifice locations are 
indicated in the tracings of Figures 2 to 7. These tracings were 
obtained from the casting of the ice shape instead of from the 
tunnel ice accretion. This was accomplished by fabricating an 
additional 15-in. casting section that was cut at the tap location 
for tracing. 

4.0 Results and Discussion 

4.1 Clean Airfoil Aerodynamics 

The NACA 23012 airfoil model was tested in the baseline 
configuration to establish the clean performance over the 
given range of Reynolds and Mach numbers. Overall, the 
clean-airfoil results followed classic airfoil behavior and 
compared favorably with existing data. Validation of the clean 
performance was established primarily by comparison of the 
present data with archival data and XFOIL results for the 
NACA 23012 airfoil section. An example of this comparison 
is shown in Figure 9 for approximately matched Reynolds and 
Mach number conditions. Broeren et al. (Ref. 5), performed a 
more recent test of a 36-in. chord NACA 23012 airfoil in the 
NASA LTPT at Re = 7.5><10 6 and M = 0.21. Abbott and von 
Doenhoff (Ref. 26) tested a 24-in. chord NACA 23012 airfoil 
also at LTPT. These data are for Re = 8.8 X 10 6 and M < 0.17 
(the exact Mach number was not given). The lift coefficient 
versus angle of attack trend shows very good agreement. The 
Broeren et al. (Ref. 5) data had a slightly higher slope for 
angles of attack preceding the stall, while the Abbott and von 
Doenhoff (Ref. 26) data had a slightly lower slope for negative 
angles of attack. The sharp stall of the clean airfoil is 
indicative of leading-edge stall. For this stall type, boundary- 
layer separation occurs near the leading edge without 
subsequent reattachment, resulting in separated flow over the 
airfoil and the significant decrease in lift (Ref. 27). The 
agreement in maximum lift coefficient was very good and well 
within the range of experimental uncertainty, potential 
variations in wall-correction methods, and manufacturing 


TABLE 3.— ICE ACCRETIONS SELECTED FOR AERODYNAMIC TESTING 


Ice 

classification 

Run 

no. 

Airspeed 

knots 

a 

deg. 

MVD 

\im 

LWC 

g/m 3 

Total 

temperature 

°F/°C 

Static 

temperature 

°F/°C 

Spray 

time 

min. 

Horn 

EG1164 

175 

5.0 

20 

0.85 

28.0/-2.2 

20.8/-6.2 

11.3 

Streamwise 1 

EG1162 

150 

2.0 

30 

0.55 

-8.0/-22.2 

-13.5/-25.3 

10.0 

Roughness 1 

EG1126 

200 

2.0 

20 

0.50 

28.0/-2.2 

18.6/-7.4 

2.0 

Spanwise Ridge 

EG1159 

150 

1.5 

20 

0.81 

20.0/-6.7 

14.8/-9.6 

15.0 

Streamwise 2 

EG1125 

200 

2.0 

15 

0.30 

4.0/-15.6 

-5.3/-20.7 

20.0 

Roughness 2 

EG1134 

200 

2.0 

40 

0.55 

4.0/-15.6 

-5.3/-20.7 

2.0 
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Figure 2. — Tracing and photograph of horn shape EG1164. 
Tracing taken from ice casting used for aerodynamic testing 
with pressure orifice locations indicated by open circles. 


Figure 3. — Tracing and photograph of streamwise shape 1 
EG1162. Tracing taken from ice casting used for 
aerodynamic testing with pressure orifice locations indicated 
by open circles. 
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Figure 4. — Tracing and photograph of roughness shape 1 
EG1126. Tracing taken from ice casting used for 
aerodynamic testing with pressure orifice locations indicated 
by open circles. 


Figure 5. — Tracing and photograph of spanwise-ridge shape 
EG1159. Tracing taken from ice casting used for 
aerodynamic testing with pressure orifice locations indicated 
by open circles. 
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Figure 6. — Tracing and photograph of streamwise shape 2 
EG1125. Tracing taken from ice casting used for 
aerodynamic testing with pressure orifice locations indicated 
by open circles. 


Figure 7. — Tracing and photograph of roughness shape 2 
EG1134. Tracing taken from ice casting used for 
aerodynamic testing with pressure orifice locations indicated 
by open circles. 
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Figure 8. — Completed installation of horn shape EG1164 
casting sections on the leading edge of the NACA 23012 
airfoil model (upper-surface view). 
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Figure 9. — Comparison of clean NACA 23012 airfoil section 
performance from the current work with results from Broeren 
et al. (Ref. 5), Abbott and von Doenhoff (Ref. 26) and XFOIL 
(Ref. 28) at closely matched Reynolds and Mach number. 
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tolerances of the various models. The present data had 
C/ max = 1.82 while the values for Broeren et al. (Ref. 5) and 
Abbott and von Doenhoff (Ref. 26) were 1.81 and 1.80, 
respectively. There was a slightly larger discrepancy in 
stalling angle of attack with a sta u = 18.1° for the present data 
and 17.6° and 18.4° for the other sources. Some of this 
disagreement could be due to the angle of attack resolution at 
which the data were acquired and flow angularity in the wind- 
tunnel. For example, the Broeren et al. (Ref. 5) data were only 
acquired in one degree increments, so it is quite possible that a 
higher stalling angle could have been measured with finer 
resolution in that test. The XFOIL results were computed for a 
Reynolds number of 9.0 X 10 6 and a Mach number of 0.20. 
XFOIL is an airfoil analysis code that couples a panel method 
flowfield solver to an integral boundary-layer formulation 
(Ref. 28). The lift coefficient trends show that the XFOIL 
results compared favorably with all of the experimental data. 
The former had a slightly higher lift-curve slope which is a 
common characteristic of XFOIL results, in the authors’ 
experience. 

The agreement in the quarter-chord pitching moment 
variation with angle of attack was very good between the 
present data and the Broeren et al. (Ref. 5) data, with the latter 
being slightly more nose down for positive angles of attack. 
The agreement with the Abbott and von Doenhoff (Ref. 26) 
pitching-moment data is not as good owing to some noted 
deficiencies in their balance arrangement (Refs. 5 and 26). The 
XFOIL calculations also reproduce the appropriate C„, versus 
a trend, matching up best with the Broeren et al. (Ref. 5) data. 
For drag coefficient, the agreement between the present data 
and Abbott and von Doenhoff (Ref. 26) is good, particularly 
for lift coefficients less than about 1.0. Both data sets have the 
characteristic reduction in C,i at C/ ~ 0.10, also present in the 
XFOIL results. The large difference in the Broeren et al. (Ref. 
5) Cd data was shown to be related to the mismatch of 
spanwise running seams for the removable leading edge on 
that model (Ref. 29). While the present model also had a 
removable leading edge, the match of the spanwise running 
seams was excellent. The Abbott and von Doenhoff (Ref. 26) 
model was solid, i.e., no seams. The XFOIL results agreed 
fairly well with the present data and Abbott and von Doenhoff 
for most of the lift range until C/ ~ 1.0. These performance 
comparisons show that valid data were obtained for the clean 
NACA 23012 airfoil section in the current work. 

The pressurization capability of the FI facility allowed for 
the independent exploration of Reynolds and Mach number 
effects on the airfoil performance. The performance results for 
the NACA 23012 exhibited classic airfoil behavior. For 
example, the lift data show that Ci„ wx increased from 1.76 at 
Re = 4.6xl0 6 to 1.84 atRe = 8.1xl0 6 to 1.88 at Re = 12.3xl0 6 
with M = 0.10. The stalling angle also increased over this 
range of Reynolds number and the pitching moment values 
near stall also became slightly more nose down. These 
Reynolds number trends agree with archival data for the 
NACA 23012 and data for other airfoils. The effect on 
maximum lift coefficient is summarized in Figure 10. Included 
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Figure 10. — Comparison of the effect of Reynolds number on 
maximum lift coefficient for the NACA 23012 airfoil from the 
current work with Broeren et al. (Ref. 5), Abbott and von 
Doenhoff (Ref. 26), NACA 0012 airfoil data from Ladson (Ref. 
30), NLF 0414 airfoil data from Addy (Ref. 9) and GLC-305 
airfoil data from Addy et al. (Ref. 6). 


in Figure 10 are NACA 23012 section data from Broeren et al. 
(Ref. 5) Abbott and von Doenhoff (Ref. 26), NACA 0012 data 
reported by Ladson (Ref. 30), NLF-0414 data from Addy 
(Ref. 9), and GLC-305 data from Addy et al. (Ref. 6). The 
NACA 0012 data were acquired in the LTPT on a 24-in. chord 
model. The NLF-0414 and GLC-305 data were also acquired 
in the LTPT using 36-in. chord models. The GLC-305 section 
was fairly thin (less than 9 percent) and not highly cambered, 
thus resulting in the lower C/ max values shown in Figure 10. 
These data indicate that large increases in Q max occur for 
Reynolds numbers less than about 6.0 X 10 6 and increasing 
Reynolds number above this results in only modest increases 
in Ci, nax . These results compare favorably with Haines’s 
(Ref. 31) detailed review of scale and Reynolds number 
effects on airfoil stalling characteristics. There is good 
agreement in the C/ max values for the three different sources of 
NACA 23012 section data in Figure 10, especially given the 
difference in models, model scales, facilities and installations. 
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Figure 11. — Comparison of the effect of Mach number on 
maximum lift coefficient for the NACA 23012 airfoil from the 
current work with Broeren et al. (Ref. 5), NACA 0012 airfoil 
data from Ladson (Ref. 30), NLF-0414 airfoil data from Addy 
(Ref. 9) and GLC-305 airfoil data from Addy et al. (Ref. 6). 


Some of the difference can also be attributed to differences in 
Mach number. 

The effects of compressibility were documented for the 
clean NACA 23012 airfoil section by varying Mach number 
from 0.10 to 0.28 at a constant Reynolds number of 12.1xl0 6 . 
Once again, classic airfoil behavior was observed in the 
characteristic increases in the lift-curve slope and pitching- 
moment slope just prior to stall for increasing Mach number. 
The maximum lift coefficient and stalling angle of attack both 
decreased with increasing Mach number, as expected. These 
results for C/ max are summarized in Figure 11 along with data 
from the other sources mentioned in connection with Figure 
10. The data for the NACA 0012 and GLC-305 airfoils show 
more significant reduction in C/ max for Mach numbers larger 
than 0.20, while the data for the NACA 23012 and NLF-0414 
sections show mixed results. The sensitivity to Mach number 
is linked to local supersonic flow near the leading edge at high 
angles of attack, and so depends, in part, on the airfoil’s 
leading-edge radius of curvature. The present data indicate a 


more linear decrease in Q max over the given Mach number 
range, while the trend in the Broeren et al. (Ref. 5) data aligns 
more closely with the other airfoils. Despite these minor 
differences, the overall Reynolds and Mach number effects 
and the comparison to historical data yield confidence in these 
measurements. 

4.2 Iced- Airfoil Aerodynamics 

The ice-casting simulations applied to the NACA 23012 
airfoil caused a large range of detrimental performance effects 
due to the differences in the ice accretion geometries and 
roughness levels. The general range in performance 
degradation is illustrated in Figure 12 for four of the six 
artificial ice shapes at Re = 15.9><1 0 6 and M = 0.20. The clean 
airfoil Ci max value was 1.85 and a sta u was 18.1° for this 
Reynolds and Mach number condition. The EG 11 62 
streamwise and EG1126 roughness shapes had a very similar 
effect on lift, drag and pitching moment. The effect of the 
EG 1 162 ice simulation was to reduce C/,„ m to 1.16 and «,,„// to 
11.9°. The EG1126 roughness shape caused a slightly larger 
lift penalty, yielding Q mar = 1.09 at a staU = 11.4°. Both of 
these ice simulations caused a significant redistribution of 
surface pressure resulting in an increase in the pitching- 
moment slope. The stall characteristics and surface pressure 
distribution indicate that there was likely some boundary-layer 
separation on the aft portions of the airfoil. The sharp drop in 
Ci near Q ma * is reminiscent of the clean leading-edge stall 
type. The plot of drag coefficient shows that the EG1162 and 
EG1126 shapes also had a similar effect, with the former 
having a slightly higher C d up to a = 8.5°, where there was a 
crossover. At higher angles of attack the EG 11 62 streamwise 
ice shape caused a lower C d , as compared to the EG 1126 
shape, which is consistent with the higher lift coefficients. 

The performance results with the streamwise (EG 11 62) and 
roughness (EG1126) shapes sharply contrast the effect of the 
horn (EG1164) and spanwise-ridge (EG1159) shapes. The size 
and location of the latter two ice shapes on the airfoil caused 
large upper-surface separation bubbles that significantly 
altered the flowfield and the subsequent performance. The 
resulting maximum lift coefficient for the EG 1 164 horn shape 
was 0.86 at a sm n = 8.8° which amounted to a 54 percent 
reduction in Q mar from the clean airfoil at this Reynolds 
and Mach number condition. The effect of the EG1159 
spanwise-ridge shape was even more severe with Q ma * = 0.52 
at a sta ii = 5.6°. The lift and pitching moment data for the airfoil 
with these ice shapes are consistent with the understanding of 
the flowfield effects relating to the large separation bubble 
formed downstream of the horn and ridge (Refs. 13 and 14). 
As discussed at length by Bragg et al., (Refs. 13 and 14) the 
separation bubble is known to have large-scale unsteady 
characteristics. In the time-averaged sense, the separation 
bubble grows rapidly in chordwise extent with increasing 
angle of attack, until the unsteady reattachment region 
approaches the airfoil trailing-edge. The time-averaged bubble 
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Figure 12. — Comparison of performance effects of the 
simulated ice configurations on the NACA 23012 airfoil at 
Re = 15.9x10® and M= 0.20. 

growth is a feature of the classic thin-airfoil stall type 
described in detail by McCullough and Gault (Ref. 27). The 
plot of drag coefficient in Figure 12 further illustrates the 
significant impact of the separation bubble for the EG 11 64 
and EG1 159 ice-shape simulations. It is interesting to note that 
C,i for the EG1159 spanwise ridge is lower than that for the 
EG1164 horn shape for angles of attack less than about -2.5°. 
This most likely occurred because at low angle of attack drag 
coefficient is affected primarily by ice accretion on the airfoil 
lower surface. A comparison of Figures 2 and 5 clearly shows 
that the horn shape had larger height and extent of ice 


accretion on the lower surface, albeit closer to the leading 
edge. This effect of lower and upper-surface ice accretion on 
drag is described by Bragg et al. (Ref. 14) and in much more 
detail by Kim (Ref. 32). For angles of attack larger than -2.5°, 
the drag coefficient for the spanwise-ridge shape was much 
larger, which was consistent with the lower C imax . The fact 
that the spanwise-ridge ice simulation resulted in the largest 
performance degradation was likely due to the location of the 
upper- surface ridge. As depicted in Figures 2 and 5, the horn 
height was larger than the upper-surface ridge, but the ridge 
was located at xlc ~ 0.06. Lee and Bragg (Ref. 33) showed that 
for the NACA 23012 airfoil, the most sensitive location for ice 
accretion (of this size) in terms of loss in maximum lift was 
downstream of the leading edge. This is also illustrated in 
Figure 12 for the streamwise shape that formed on the airfoil 
leading edge, as the performance penalty was not nearly as 
severe as for the spanwise ridge. 

The effect of leading-edge ice accretion on the NACA 
23012 airfoil performance is further illustrated in Figure 13. 
These data are for the EG1125 and EG1134 streamwise and 
roughness ice shapes compared to the streamwise and 
roughness shapes plotted in Figure 12. Of these, it is 
interesting to note that the lowest iced-airfoil Q max occurred 
with the EG1126 roughness shape. The drag penalty was also 
largest for this shape for angles of attack greater than 9°. Note 
that there is a scale change between Figures 12 and 13 in the 
C,i versus a plot to facilitate the comparisons. The fact that the 
two streamwise ice simulations (EG 11 62 and EG 1125) had a 
nearly identical effect on Q mal and C,j (for a > 9°) as the 
EG1126 roughness shape reveals an important characteristic 
about leading-edge ice accretion on the NACA 23012 section. 
A comparison of the ice accretion characteristics in Figures 3, 
4 and 6 and icing environment in Table 2, shows a wide range 
of icing conditions and resulting sizes and shapes. The 
obvious implication is that the basic ice height or thickness for 
these types of accretions is not a strong driver of the stall 
performance behavior. A similar effect has been observed in 
previous testing of a subscale NACA 23012 airfoil with larger 
geometric quarter-round ice simulations located at the leading 
edge (i.e., xlc =ylc = 0). Lee and Bragg (Ref. 33) showed that 
the iced airfoil Ci f , nax varied from about 1.01 to 0.97 as the ice 
height (or thickness) varied from k/c = 0.006 to 0.014. A 
similar Q„ MA range was observed in Figure 13; 1.09 for the 
EG1126 simulation to 1.16 for the EG1162 simulation. 
Despite the similarity in Ci max values, there is a difference in 
the stalling characteristics with the EG 1162 and EG 1126 
simulations having a sharper, leading-edge type stall versus 
the EG1125 simulation that resulted in a more gradual trailing- 
edge type stall. This implies that significant flowfield 
differences remain among the various artificial ice shapes 
despite the alignment of maximum lift coefficient. These 
effects are addressed further in Section 4.4. 

The data for the EG1134 simulation illustrate the smallest 
maximum lift penalty measured for the given set of artificial 
ice shapes. The C !max of 1.28 was about 10 percent higher than 
for the other three ice simulations in Figure 13. The drag 
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Figure 13. — Comparison of performance effects of the 
roughness and streamwise ice simulations on the NACA 
23012 airfoil at Re = 15.9x10® and M = 0.20. 


coefficient on the interval —1.5 < a < 4.7° was slightly higher 
than, but most similar to the drag coefficient for the other 
roughness simulation, EG1126. It is interesting that the lift 
performance of the NACA 23012 section with the EG1134 
simulation was very similar to the 80-grit sandpaper 
(k/c = 0.00023) reported by Broeren et al. (Ref. 5) on the 36- 
in. chord model of the same airfoil. For that case, the Q max 
was about 1.30 at Re = 7.5><1 0 6 and M= 0.21. The minimum 
G c i for the airfoil with the 80-grit sandpaper was 0.010, about 
10 percent higher than for the EG1134 roughness simulation 


in the present data. The comparison of the EG1134 
performance effects to the 80-grit sandpaper on the 36-in. 
chord model is useful because the latter is roughly equivalent 
to 40-grit sandpaper on a full-scale wing, which is often used 
in aircraft certification flight testing applications (Ref. 34). 

The trends in the drag coefficient variation with angle of 
attack for the iced-airfoil configurations in Figure 13 may be 
interpreted in terms of the icing environment from which these 
simulations were obtained. In the case of the two streamwise 
ice simulations, for example, the drag coefficient for the 
EG1162 configuration was higher than for the EG1125 
configuration in the range of -4° to 8°. This is somewhat 
surprising given that the former ice shape was more conformal 
to the airfoil leading edge than the latter. Flowever, the former 
was accreted with a larger MVD drop size and therefore the 
downstream roughness covered a larger surface extent. This 
may be at least partially responsible for the larger drag over 
this range. It was noted above that the drag coefficient for the 
EG1134 roughness configuration was slightly higher than for 
the EG1126 roughness configuration over the interval —1.5 < 
a < 4.7°. This may have occurred because the rime roughness 
EG1 134 was accreted with a larger MVD drop size resulting in 
a larger surface extent of ice roughness. For angles of attack 
larger than 4.7°, it was likely that a larger extent of trailing- 
edge separation developed for the EG1126 configuration, thus 
the drag coefficient grew larger than for the EG1134 
configuration corresponding to the lower stalling angle. 

4.3 Iced-Airfoil Reynolds and Mach Number 
Effects 

The effect of Reynolds number variation at constant Mach 
number was explored for all of the ice-shape configurations. 
The pitch-polar data for the EG 1125 streamwise-ice 
simulation are shown and discussed in Broeren et al. (Ref. 35). 
The lift and pitching moment data were virtually unchanged 
despite a nearly three-fold increase in Reynolds number from 
4.6 X 10 6 to 12.0 X 10 6 . There was a small increase in maximum 
lift coefficient from 1.10 at the lowest Reynolds number to 
1.13 at the highest. The pitching moment data show that there 
was a small 0.5° increase in the stall break as the Reynolds 
number was increased from 4.6 X 10 6 to 8.2 X 10 6 . The drag data 
show slightly more dependence with the drag coefficient 
generally decreasing with increasing Reynolds number. The 
reason for this trend was not investigated in detail, but this 
effect on drag has been observed in other iced-airfoil studies 
(Ref. 5 and 6). 

The Reynolds and Mach number effects data for the 
EG1125 configuration was representative of all of the ice 
simulations tested in this study. This was particularly true for 
the lift and pitching moment variations with Reynolds number, 
even for the small roughness shapes. Changes in Reynolds 
number over this range did not substantially change the 
character of the stall. That is, iced configurations with an 
abrupt, leading-edge stall character maintained this over the 
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Figure 14. — Effect of Reynolds number on maximum lift 
coefficient for the clean NACA 23012 airfoil and with the 
simulated ice shapes. 

test Reynolds number range. The effect on maximum lift 
coefficient for all of the simulations tested is summarized in 
Figure 14. Also included in the figure are the data for 
M = 0.20. The non-effect of Reynolds number is consistent 
with the previous data compiled by Bragg et al. (Ref. 14). The 
authors argue that the presence of ice contamination on the 
airfoil governs the boundary-layer behavior such that the 
influence of Reynolds number on such processes is reduced to 
lower order effects unlike in the clean case. The present data 
continue to support this conclusion. 

The data in Figure 14 were combined with data compiled 
from other sources in Figure 15. The purpose of this figure is 
to further illustrate the difference in clean versus iced-airfoil 
Reynolds number effects on maximum lift coefficient, which 
is an important aircraft safety-related parameter. Incorporated 
in Figure 15 are data for four different airfoils with numerous 


-• Present Data, NACA 23012, Clean, M= 0. 10 

^ Present Data, NACA 23012, Clean, M= 0.20 

-m Broeren et al., NACA23012, Clean, M = 0. 12 

♦ Ladson, NACA 0012, Clean, M= 0.20 

H Addy, NLF-0414, Clean, M= 0.12 

— «t Present Data, NACA 23012, L.E. Ice, M=0.10 
-A Present Data, NACA 23012, L.E. Ice, M= 0.20 


-a Broeren et al., NACA 23012, Roughness, M = 0.12 

-0 Broeren et al, NACA 23012, Intercycle, M= 0.12 

* & Addy, NLF-0414, L.E. Ice, M= 0. 12 



Reynolds Number, Re (xlO 6 ) 


Figure 15. — Comparison of the effect of Reynolds number on 
maximum lift coefficient for the NACA 23012 airfoil from the 
current work with Broeren et al. (Ref. 5), NACA 0012 airfoil 
data from Ladson (Ref. 30), NLF 0414 airfoil data from Addy 
(Ref. 9) and GLC-305 airfoil data from Addy et al. (Ref. 6). 

different kinds of simulated leading-edge ice accretion. These 
data show that for some iced-airfoil cases, most notably the 
sandpaper roughness due to Broeren et al. (Ref. 5), there can 
be measurable Reynolds number effects on C/ max for Re < 
4.0 X 10 6 . The other data shown in Figure 15 that were taken 
from Broeren et al. (Ref. 5) are for “intercycle” ice shapes. 
These ice shapes, being larger, leading-edge type ice, did not 
result in significant Q m „ decreases for Re < 4.0 X 10 6 . There is 
at least one other study (Ref. 36) (not plotted here) that further 
shows negligible variation in iced-airfoil Ci max from 
Re =2.0 X 10 6 to 10.5 X 10 6 . Thus, increases in Reynolds number 
above 2.0 X 10 6 tend to have little or no effect on C/ max for 
airfoils with large, leading-edge ice accretion. This is an 
important result for subscale simulation at low-Reynolds 
number. 
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Mach Number, M 

Figure 16. — Effect of Mach number on maximum lift coefficient 
for the clean NACA 23012 airfoil and with the simulated ice 
shapes. 

Broeren et al. (Ref. 35) also show the effect of Mach 
number variation at constant Reynolds number in detail for the 
same EG 1125 streamwise-ice simulation. The data show that 
the slopes of C/ and C,„ with respect to angle of attack were 
increased prior to stall as the Mach number was increased 
from 0.10 to 0.28. The maximum lift coefficient was reduced 
from 1.13 to 1.11 and stalling angle was reduced from 11.9° to 
10.9° over this range. The minimum drag coefficient increased 
by a larger amount for Mach number increasing from 0.10 to 
0.20, with little change for Mach number increasing from 
0.20 to 0.28. In contrast, drag coefficients at higher angles of 
attack were influenced more strongly for increasing Mach 
number from 0.20 to 0.28. Mach number performance sweeps 
for the other iced-airfoil configurations yielded similar results. 
The effect on maximum lift coefficient is summarized in 
Figure 16. The horn shape (EG1164) and the spanwise-ridge 
shape (EG1159) showed the most significant decrease in C/,„ oi 
over this range, with the smaller ice shapes showing much less 
dependence of Ci max on M. This trend involving smaller ice 
shapes and Mach number effects on C tmax has been observed 
in other studies. 

Plotted in Figure 17 are data compiled from a number of 
iced-airfoil studies. For example, Broeren et al. (Ref. 5) report 
data for two sandpaper roughness configurations and four 


-• Present Data, NACA 23012, Clean, Re = 12.1X10 6 

Broeren et at., NACA 23012, Clean, Re = 10.5xl0 6 

-■ Ladson, NACA 0012, Clean, Re = 9.0xl0 6 

-♦ Ladson, NACA 0012, Clean, Re = 6.0xl0 6 

-* Addy, NLF-0414, Clean, Re = 6.4xl0 6 

G - Present Data, NACA 23012, L.E. Ice, Re = 12.2xl0 6 

-a Broeren et al., NACA 23012, Roughness, Re = 7.5xl0 6 

-0 Broeren et al., NACA 23012, Intercycle, Re = 7.5xl0 6 

-® Addy, NLF-0414, L.E. Ice, Re = 6.4xl0 6 

« Addy et al., GLC-305, L.E. Ice, Re= 10.5xlQ < ’ 



Figure 17. — Comparison of the effect of Mach number on 
maximum lift coefficient for the NACA 23012 airfoil from the 
current work with Broeren et al. (Ref. 5), NACA 0012 airfoil 
data from Ladson (Ref. 30), NLF-0414 airfoil data from Addy 
(Ref. 9) and GLC-305 airfoil data from Addy et al (Ref. 6). 


intercycle ice configuration on the NACA 23012 section. For 
the two sandpaper roughnesses, there was a very slight 
decrease in C/ max for one case and an increase in C/, max for the 
other case over a Mach number range of 0.12 to 0.28. In 
contrast, the much larger, ridge-type intercycle ice 
configurations resulted in Ci max variations with Mach number 
similar to that shown in Figure 16 for the EG 1164 and 
EG1159 simulations. The data for the leading-edge, iced- 
airfoil configurations on the NLF-0414 and GLC-305 airfoils 
from Addy et al., (Refs. 6, 7 and 9) also exhibit similar Mach 
number trends. The decrease in Ci„ lCLX with increasing Mach 
number for larger ice shapes appears to be related to an 
increase in the size of the separated flow region and is 
described in more detail by Bragg et al. (Ref. 14) and Broeren 
et al. (Ref. 36). 
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4.4 Discussion of Iced- Airfoil Stalling 

Characteristics 

The stalling characteristics of the iced-airfoil configurations 
tested in this study are of key interest, not only as they relate 
to the safety of flight operations, but also as they relate to the 
subscale component of this work described by Bragg et al. 
(Ref. 12). The lack of significant Reynolds and Mach number 
effects on maximum lift coefficient implies that simulation of 
iced performance at small scale may be quite successful. In a 
given subscale simulation, it is also important that the stalling 
characteristics are maintained in addition to any one parameter 
such as C'i max or a sta u. The development of airfoil stall 
classifications relates back to early research of B. Melvill 
Jones (Ref. 38), among others. McCullough and Gault 
(Ref. 27) built upon the foundations of this early work and 
conducted a systematic series of airfoil section tests to 
establish the commonly accepted stall-type definitions in use 
today. Trailing- edge stall is preceded by movement of the 
turbulent boundary-layer separation point forward from the 
trailing edge with increasing angle of attack. Leading-edge 
stall has abrupt flow separation near the leading edge 
generally without subsequent reattachment. The “abrupt” 
separation usually results from the “bursting” of a small 
laminar separation bubble that results in a sharp decrease in 
lift. Thin-airfoil stall is preceded by flow separation at the 
leading edge with reattachment (laminar separation bubble) at 
a point that moves progressively downstream with increasing 
angle of attack. Stall type is a function of several variables 
such as Reynolds number, surface roughness or free-stream 
turbulence. Therefore, any particular airfoil may exhibit a 
combination of stall types, or its stall type may change over 
various flow regimes or conditions. 

Over the range of Reynolds and Mach numbers in this study 
the clean NACA 23012 airfoil clearly exhibits the leading- 
edge stall type. The abrupt flow separation at stall was evident 
from the flow visualization and pressure distributions and is 
manifest in the sharp drop in lift coefficient. Both Haines 
(Ref. 31) and Tani (Ref. 39) conducted extensive reviews of 
airfoil stall research and noted that the abrupt flow separation 
associated with leading edge stall can result from either the 
bursting of a small laminar separation bubble (as reported by 
McCullough and Gault (Ref. 27)) or from re-separation of the 
turbulent boundary-layer downstream of the bubble 
reattachment location. Haines makes the case that the latter 
separation mechanism is probably more likely for moderate to 
high-Reynolds number flows. So the bursting of the laminar 
separation bubble is not a necessary prerequisite for leading- 
edge stall. In either case, there is no question about the basic 
stall type of the clean NACA 23012 section. 

The application of the simulated ice shapes to the NACA 
23012 airfoil section alters the stall type in some cases. It is 
clear that the main stall mechanism for the EG 1164 horn 
shape and EG1159 spanwise-ridge shape configurations was 
the growth (with angle of attack) of a large upper-surface 
separation bubble. Thus, the airfoil with these ice shapes 


would be classified as having the thin-airfoil stall type. This 
conclusion is based on lift performance data in Figure 12 and 
corroborated by the corresponding pressure distributions and 
flow visualization (not shown). These flowfield physics are 
fairly well documented (e.g., see Bragg et al. (Ref. 13 and 
14)). However, the stalling characteristics of the streamwise 
and roughness ice simulations on the NACA 23012 airfoil are 
not as well understood and deserve a closer look. 

What is of primary interest here is that for three of the four 
roughness and streamwise ice simulations tested on the NACA 
23012 airfoil, the presence of the simulated ice did not appear 
to significantly alter the leading-edge stall type. This 
observation is based upon the abrupt loss in lift at stall for the 
EG1162, EG1126 and EG1134 simulations (cf. Fig. 13). Only 
the EG1125 simulation resulted in a gradual loss of lift post 
stall. This comparison is illustrated in more detail in the 
pressure distributions where the EG1134 configuration is used 
as the example for leading-edge stall. The data plotted in 
Figure 18 for the iced-airfoil case at a = 13.5° and 13.8° show 
the profound redistribution of surface pressure as the stall 
occurs with only a 0.3° increase in angle of attack. The nearly 
constant pressure region aft of about x/c = 0.10 results from 
the massive separation on the upper surface indicative of the 
leading-edge stall type. Also shown for comparison is the 
clean pressure distribution for a = 13.5° At this angle of attack 
the clean and iced pressures were very similar. There was 
some divergence of the pressure near the trailing edge in the 
latter case indicating that there was likely some extent of 
turbulent boundary-layer separation. This separation combined 
with the effects of the ice shape probably contributed to the 
reduced suction pressures on the forward portion of the airfoil 
relative to the clean configuration. Technically, the iced-airfoil 
stall type may be classified as a combination of leading-edge 
and trailing-edge stall due to the presence of the turbulent 
boundary-layer separation. However, the abrupt redistribution 
of surface pressure (at a = 13.8°) is a clear hallmark of 
leading-edge stall. 

This behavior is contrasted with the data in Figure 19 for 
the EG1125 streamwise-ice simulation. The iced-airfoil 
pressure distributions correspond to maximum lift (a = 11.9°) 
and 1.5° post stall (a = 13.3°). In contrast to the EG1134 case, 
the pressure distributions are comparable corresponding to the 
gradual reduction in lift coefficient shown in Figure 13. There 
was more of a divergence of C p on the aft section of the airfoil 
at the higher angle of attack indicating a larger extent of 
turbulent boundary-layer separation. However, the salient 
feature for both angles of attack is the near-constant pressure 
region between x/c = 0.015 and 0.05. This probably indicates 
the presence of a separation bubble induced by the simulated 
ice shape. The bubble significantly altered the pressure 
distribution from the clean case, such that the tendency to 
leading-edge stall was altered. In fact, the EG1125 
configuration could be classified as having trailing-edge stall, 
since the reduction in lift occurred due to the turbulent 
boundary-layer separation moving forward as angle of attack 
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Figure 18. — Chordwise pressure distribution comparison for 
angles of attack near stall for the EG1134 roughness 
simulation on the NACA 23012 airfoil at Re = 15.9x10 6 and 
M= 0.20. 



Figure 19. — Chordwise pressure distribution comparison for 
angles of attack near stall for the EG1 125 streamwise shape 


on the NACA 23012 airfoil at Re 
M= 0.20. 


15.9x10° and 


increased. While there was a small separation bubble present 
due to the ice shape, this bubble did not grow significantly in 
chordwise extent leading up to the stall. This partly explains 
the classification of this ice shape as streamwise ice, according 
to the descriptions established by Bragg et al. (Refs. 13 and 
14). This behavior is contrasted with the EG1 164 and EG1 159 
configurations that had the thin airfoil stall type. The stall 
mechanism for the EG1125 configuration is very different 
from the EG1164 and EG1159 configurations even though 
they shared a more gradual decrease in lift at stall. This is 
consistent with the characteristics of thin-airfoil and trailing- 
edge stall as defined by McCullough and Gault (Ref. 27). 

The fact that the leading-edge stall characteristics of the 
NACA 23012 airfoil were not significantly altered by the 
presence of two roughness and one streamwise-ice shapes is 
important for subscale simulations and computational 
modeling. As mentioned in this discussion, leading-edge stall 
is defined as abrupt flow separation near the leading edge of 
the airfoil without subsequent reattachment. With the ice 
simulations in question, there were likely small areas of 
boundary-layer separation present on the scale of the ice 
roughness. Research reported by Kerho and Bragg (Ref. 40) 
showed that leading-edge roughness can lead to a 
“transitional” boundary layer instead of acting as a “trip” with 
energetic transition to turbulence. The transition process due 
to ice-type roughness develops over a large part of the airfoil 
chord, up to 30 to 50 percent in some cases. The present data 
indicate that the “transitional” nature of the boundary layer is 
still susceptible to abrupt separation for the three ice 
simulations. In contrast, the other ice shapes being larger, 
caused larger-scale separation and alteration of the transition 
process as well as the pressure distribution, thus leading to an 
alteration of the stall type. 

Another important feature of the roughness and streamwise 
ice shapes was that the effect on maximum lift was similar in 
magnitude despite the large variation in geometry. As 
discussed in connection with Figure 13, the range in Q max was 
from 1.09 to 1.28. While this variation is not insignificant, it is 
small compared to the decrease from the clean value of 1.85. 
These data were compared to other results for simulated 
leading-edge ice contamination on the NACA 23012 airfoil 
shown in Figure 20. There appears to be some dependence of 
Cl, max on ice (roughness) size located at the leading edge up to 
k/c ~ 0.0015. For sizes larger than this, there is no clear trend 
among the data gathered from various sources. The selected 
data are for a large range of ice-shape geometries and 
simulation methods. The roughness cases covered the leading 
edge, extending downstream several percent chord. The 
quarter-round data were simple geometric shapes with the 
forward face located at the leading edge. The intercycle 322 
shape had characteristics similar to the EG 11 62 ice shape in 
the current study. The reduced sensitivity to these leading- 
edge ice simulations may be due in part to the favorable 
pressure gradient in this region. Since there is no clear 
correlation to the height of the simulated ice, the geometry, 
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Figure 20. — Comparison of Ci, max dependence on ice 
(roughness) height for various simulations all applied to the 
leading-edge region of the NACA 23012 airfoil. In addition to 
the present data, the 80- and 40-grit sandpaper and 
intercycle 322 simulations are from Broeren et al. (Ref. 5). 
The NACA standard roughness datum is from Abbott & von 
Doenhoff (Ref. 26). The quarter-round data are from Lee 
and Bragg (Ref. 33). 


surface extent and distribution of the roughness play an 
important role in determining maximum lift. This type of 
result, while perhaps obvious in light of the present data, 
supports the notions of Bragg et al. (Refs. 13 and 14) in their 
development of the ice-shape classifications. This leads to an 
ironic conclusion for subscale simulation. The reduced 
sensitivity may allow for lower-fidelity simulations to yield 
acceptable values of maximum lift coefficient. The 
disadvantage is that the physics of the flowfield may not be 
properly modeled in the subscale case. Therefore, it may 
actually be more difficult to develop proper subscale 


simulations due to the reduced performance sensitivity of the 
roughness and streamwise ice shapes. 

5.0 Conclusions 

The purpose of this article is to report the results of a study 
establishing a high-fidelity, full-scale, iced-airfoil 
aerodynamic performance database. This research was 
conducted as a part of a larger program with the goal of 
developing subscale simulation methods for ice accretion. 
Airfoil performance testing was carried out at the ONERA F 1 
pressurized wind tunnel using a 72-in. (1828. 8-mm) chord 
NACA 23012 airfoil over a Reynolds number range of 
4.5xl0 6 to 16.0xl0 6 and a Mach number range of 0.10 to 0.28. 
A total of six ice-shape simulations were tested: one horn 
shape, one spanwise-ridge shape, two streamwise shapes and 
two roughness shapes. These artificial ice shapes were high- 
fidelity castings made from molds obtained during earlier ice 
accretion testing at the NASA Glenn Icing Research Tunnel. 

The artificial ice shapes had a large detrimental effect on 
the performance of the NACA 23012 airfoil. The spanwise- 
ridge shape caused the largest reduction in maximum lift with 
a value of 0.52 compared to the clean value of 1.85 at Re = 
15.9x1 0 6 and M = 0.20. The stalling angle was also reduced 
from 18.1° down to 5.6°. There was a correspondingly large 
increase in the angle of attack dependence of the pitching- 
moment coefficient. The minimum drag coefficient was 
increased by more than a factor of four from 0.0052 to 0.0224. 
The performance of the airfoil with the horn shape was better 
than with the spanwise ridge, but still severely compromised 
from the clean configuration. For the same Reynolds and 
Mach numbers, the maximum lift coefficient was 0.86 at 8.8° 
angle of attack. The minimum drag was increased by a factor 
of three from the clean value. The two roughness and 
streamwise-ice simulations had a similar range of performance 
effects on the airfoil. The range of maximum lift coefficient 
was 1.09 to 1.28 and the range of minimum drag coefficient 
was 0.0082 to 0.0106 at Re = 15.9*10 6 and M = 0.20. The 
range of performance effects was small relative to the large 
differences in the size and geometry of the roughness and 
streamwise-ice shapes tested. 

While the ice simulations negatively affected the stall 
performance of the NACA 23012 airfoil, half of the 
configurations tested resulted in a combination leading-edge 
and trailing-edge stall type. The abrupt loss of lift at stall 
associated with the leading-edge stall persisted for both the 
roughness and one streamwise-ice configurations, despite 
significant trailing-edge separation. This is an important detail 
for accurate subscale simulation because it means that while 
the ice shape was large enough to have a significant 
detrimental effect on the performance, the boundary layer in 
the leading-edge region was not significantly altered so as to 
prevent separation at the leading edge without subsequent 
reattachment. This also has important implications for the 
development of flight simulation models for iced-aircraft 
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upsets. Proper modeling of iced-aircraft stall dynamics is 
important for pilot training in recovery techniques using a 
flight simulator. For the remaining streamwise-ice simulation 
the stall type was changed to a trailing-edge stall with gradual 
lift loss at stall. For the horn and spanwise-ridge shapes the 
stall was characteristic of thin-airfoil stall since there was a 
large upper-surface separation bubble that grew with 
increasing angle of attack precipitating the stall at such low 
values of lift coefficient and incidence. 

The pressurization capability of the ONERA FI Facility 
was used to document the effects of Reynolds and Mach 
number on the clean and iced performance. The clean airfoil 
maximum lift coefficient increased from 1.76 to 1.88 as the 
Reynolds number was increased from 4.6><10 6 to 12.3 x 10 6 at 
constant M= 0.10. Increasing Mach number from 0.10 to 0.28 
at a constant Re = 12.1><10 6 had the effect of reducing the 
maximum lift coefficient from 1.88 to 1.78. The trends in lift 
and moment slope versus angle of attack and drag coefficient 
were consistent with classic airfoil behavior. For the airfoil 
with the ice simulations, there was virtually no measurable 
change in maximum lift coefficient over the entire Reynolds 
number range tested. Changes in Mach number had minor 
effects on maximum lift for the horn and spanwise-ridge 
shapes, but virtually no effects for the others. This lack of 
significant Reynolds and Mach number effects implies that 
subscale simulation of the iced-airfoil performance may be 
quite successful, despite the challenges associated with the 
roughness and streamwise-ice configurations’ stalling 
mechanism. 
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